Exactly massless quarks on the lattice. 
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Abstract 

It is suggested that the fermion determinant for a vector-like gauge theory with strictly 
massless quarks can be represented on the lattice as det •Mp-, where V = X(X^ X)~ 1 / 2 and 
X is the Wilson-Dirac lattice operator with a negative mass term. There is no undesired 
doubling and no need for any fine tuning. Several other appealing features of the formula 
are pointed out. 
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From the start, it was awkward to ensure masslessness of quarks in lattice QCD. 
Without masslessness there are no chiral symmetries and no massless pions. Massless 
QCD is a beautiful parameter free theory in the continuum, just like pure YM is. However, 
while pure YM on the lattice is elegant, the addition of quarks makes the lattice theory 
less appealing. The overlap [1, 2, 3] provided a solution to this state of affairs; however, 
the formalism is a bit complicated, and its elegance becomes a matter of personal taste. 
More importantly, the computational cost seems to be daunting. Domain walls [4, 5, 
6] provide an approximation to the overlap, but they sacrifice strict masslessness, so we 
are back to the more murky situation of standard formulations. Staggered fermions do 
preserve masslessness, but it is a nuisance to have so many of them, and the connection 
to continuum operators is messy. My purpose here is to present a simple formula for 
the fermion determinant which I argue is equivalent to the overlap but is much more 
attractive in appearance. Only future work can tell with any confidence whether real 
practical progress has been made, but there are some very attractive features that I wish 
to point out. 

The basic observation stems from recent work [7] on the overlap in odd dimensions. 
For d = 2k + 1 the overlap for Dirac fermions can be written as the determinant of a 
finite matrix of fixed shape. This simpler formula descends by dimensional reduction 
from d + 1 dimensions where it represented Weyl fermions in a less explicit formulation. 
Dimensionally reducing again, a simple formula for Dirac fermions in 2k dimensions is 
obtained. To make this letter relatively self-contained, I present a more direct derivation, 
involving no dimensional reductions. 

In [1] the overlap for a vector-like theory is constructed as follows: The chiral deter- 
minant is replaced at the regulated level by the overlap of two many-body states. These 
are the ground states of two bilinear Hamiltonians, 

= otfTV (1) 
with all indices suppressed. The matrices are obtained from 

*(»•>= (V _/_ m ) e> 

with H + = H(oo), H~ = H(—m ) and < m < 2. The infinite argument for H + can be 
replaced by any finite positive number, but the equations are somewhat simpler with our 
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choice [8, 9, 10]. The matrices C and B are given below: 

^ 2k 

(C) xai , yf 3j = - ^2vf[5y,x + ^U^x)) l3 - 5x, V +{i(Ul(y))ij], 
^ 2k 

{B) xol i,yf3j = -5 a p ypSxySjj ~ 5y !X + fx(U^(x))ij ~ 8 X ,„ + £ ( E/J ) ij 



x,y are sites on the lattice, a, (3 are Weyl spinor indices and i,j are color indices. The a M 
are Euclidean Weyl matrices in 2k dimensions. The overlap O, is given by 



= \<v+\v_>\ 2 , H ± \v±>=E± in \v±> . (4) 

An equivalent representation is also given in [1]: 

O = | < v+\v- > || < t + \t- > |, H ± \t± >= E± &x \t± > . (5) 

E^. , s denote minimal (maximal) energies which define the associated states, assuming 

min(max) V / o i o 

no degeneracies. 

I now replace the two sets of fermion operators, one for each factor in equation (5), 
by a single set of double their size: 



(6) 



Clearly, 

0=\<V+\V->\, Hf\V ± >= E±JV ± > . (7) 

But, using [1] again, we immediately can write down the overlap as 

oHdet 1 + r -r< g -> |, r 2W = (l »V e (iZ) s * (8) 



-1 J ' ' ' ~ yTfp 

where we assume det H ^ (this is equivalent to the absence of degeneracies mentioned 
after equation (5) above). The matrix 

V = T 2k+1 e(H-) (T 2 2k+1 = 1, e(H-) 2 = 1) (9) 

is unitary while both factors are hermitian. Actually, the absolute sign in eq. (8) is not 
needed: 

Ro Co 

det(l + V) = JJ(l + A r ) JJ|1 + A C | 2 > 0. (10) 



Here, A r = A* and A c 7^ A*. Equation (8), without the absolute value signs, is the main 
result of this letter. Noting that 

r 2k+1 H(m) = X(m) = ( B _ + c ™ B C +m ^ (11) 

and that 

(5 ?)-*• < i2 » 

where the 7^ are Euclidean Dirac matrices, we obtain the result announced in the abstract. 
X(m) is the familiar Wilson Dirac operator, and we denote X(— mo) by X. 

The inverse of the matrix whose determinant we have to take is easily seen to have 
only a single massless Dirac particle pole, so there is no unnecessary doubling [7] . In some 
sense our formula is a realization of an idea of Rebbi's [11], only some of the ingredients 
are different and there is a body of overlap work to rely on [1, 2, 3, 12]. 

One of the most salient properties of strictly massless fermions are the exact and 
robust zeros in instanton backgrounds. From the equivalence to the overlap we know 
that V should have —1 as an eigenvalue with degeneracy ^\tre(H~)\ (the overlap lattice 
topological charge is given by n top = |tre(i/ - )) and that this property is robust under 
small variations of the background. It is amusing to see this directly in the particular case 
of a single instanton or anti-instanton, where |rz t0 p | = 1- The following identities are easily 
established: 

det(l + V) = det(l + V' 1 ) = det f 1 + ^ = det(V) det(l + V). (13) 

det V 

If the number of sites is O, the dimension of V is N = R + 2C = 2 k djiQ, where cIr is the 
range of group indices. Clearly, N is even. 

det(V) = det(r 2fc+1 )det(e(ir)) = (-)f (-)f p = -1 (14) 

Inserting det(V) = —1 into equation (13) gives det(l + V) = 0, the expected result. 

Once we are willing to interpret as the Dirac operator, many applications in four 
dimensions suggest themselves. In the continuum, the spectral density of the massless 
Dirac operator per unit four- volume is supposed to concentrate at the origin to reproduce 
the chiral symmetry order parameter [13]. Attempts to show this on the lattice in the past 
were hindered by the eigenvalues not falling on a line in the complex plane, and it was 
difficult, in finite volumes, to identify and remove the effects of global topological charge, 
which is something one should do. Now it would be easy: all the eigenvalues of V reside 
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on the unit circle. Those strictly at —1 reflect global topology. Those close to —1, but not 
exactly at —1, should build up the condensate < ipip >. 

A related application is the measurement of f n using a finite-size soft-pion theorem 
as in [14]. Since chiral symmetry is exact we do not have to deal with the more complex 
situation created by the presence of an explicit symmetry breaking parameter. Anyhow, 
the interpretation of the transition obtained by tuning the mass parameter in ordinary 
Wilson-Dirac formulations, while away from the continuum limit, is likely different from a 
full restoration of chiral symmetries as first pointed out by Aoki [15]. Therefore, working 
with ordinary Wilson fermions and obtaining f n via finite size theorems is promising to be 
messy. 

Since the eigenvalue distribution of V has the potential of behaving smoothly on 
the lattice there is a chance to make some progress on the painful problem of QCD at 
finite chemical potential [16], at least in the rather interesting situation of strictly massless 
quarks. 

In addition to the overlap, several proposals to regulate chiral gauge theories main- 
tain exact gauge invariance of the absolute value of the chiral determinant [17,18]. It is 
typically proposed [17] to take a square root of |det(X)| with a mass parameter finely 
tuned to Wilson- Aoki criticality. Clearly, this is not very appealing, since det(X) changes 
sign frequently in that vicinity (it is real) and the absolute value will give trouble when 
variations with respect to the background fields are taken, for example for the purpose of 



deriving a Ward identity. However, using y det ^4^- avoids this problem and is identical 
to the absolute value of the chiral determinant in the overlap; as a bonus the exact zeros 
induced by nontrivial global topology are also maintained. Maybe there is a possibility to 
unify in this way several different approaches to the problem of regularizing chiral gauge 
theories. 

In recent work [19], it was shown that the method of [1] for using the overlap to 
count instantons gives results in perfect agreement with very different methods [20] , thus 
providing much support for both results. By varying mo in the interval (0, 2) a flow of 
^-eigenvalues on the unit circle would be induced. The eigenvalue motions on the unit 
circle have the potential of providing a pictorial insight via a scale-dependent (the scale is 
~ ^) fermionic probe of the gauge background. 

Much of what I said needs further study. Some of the needed future work is analytical: 
With the new formulae, perturbation theory may be more compact and more Feynman- 
like. Of course, one needs to check whether the natural guess j^y is indeed an appropriate 
definition of the fermion propagator, and how it is related to the continuum. Numerical 
work in two dimensions could check the flow patterns and see whether they are interesting. 
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Also, comparisons between the pure overlap and its domain wall approximations [21] might 
be facilitated by the new version of the formula. And, most importantly, it is an attractive, 
if difficult, challenge to develop an efficient numerical technique to include det in 
dynamical simulations in four dimensions, or even to compute just the inverse j^y- While 
the problem looks hard, it does seem somewhat easier than when looking at the older form 
of the overlap in equation (4). 

This work was supported in part by the DOE under grant # DE-FG05-96ER40559. 



References. 

[1] R. Narayanan, H. Neuberger, Nucl. Phys. B 443 (1995) 305. 
[2] R. Narayanan, H. Neuberger, Phys. Lett. B 302 (1993) 62. 

[3] S. Randjbar-Daemi and J. Strathdee, Phys. Lett. B348 (1995) 543; Nucl. Phys. B443 
(1995) 386; Phys. Lett. B402 (1997) 134; Nucl. Phys. B461 (1996) 305; Nucl. Phys. 
B466 (1996) 335. 

[4] D. B. Kaplan, Phys. Lett. B288 (1992) 342. 

[5] T. Blum, A. Soni, Phys. Rev. D56 (1997) 174. 

[6] V. Furman, Y. Shamir, Nucl. Phys. B439 (1995) 54. 

[7] Y. Kikukawa, H. Neuberger, |hep-lat/9707UT6[ 

[8] D. Boyanovsky, E. Dagotto, E. Fradkin, Nucl. Phys. B 285 (1987) 340; Y. Shamir, 

Nucl. Phys. B406 (1993) 90. 
[9] R. Narayanan, H. Neuberger, Phys. Lett. B 393 (1997) 360; Y. Kikukawa, R. 
Narayanan, H. Neuberger, Phys. Lett. B 399 (1997) 105. 
[10] Y. Kikukawa, R. Narayanan, H. Neuberger, |hep-lat/9705U06| (1997). 
[11] C. Rebbi, Phys. Lett. 186B (1987) 200. 
[12] A. Yamada, |hep-lat/9705040| . 

[13] T. Banks, A. Casher, Nucl. Phys. B169 (1980) 103. 

[14] U. M. Heller, H. Neuberger, Phys. Lett. 207B (1988) 189. 

[15] S. Aoki, Phys. Rev. D30 (1984) 2653. 

[16] I. M. Barbour, S. E. Morrison, E. G. Klepfish, J. B. Kogut, M-P Lombardo, |Eep 
| lat/ 9705038 . 



[17] S. Delia Pietra, V. Delia Pietra, L. Alvarez-Gaume, Commun. Math. Phys. 109 
(1987) 691; K. Okuyama, H. Suzuki, |hep4at/9706163| ; G. T. Bodwin, Phys. Rev. 



D54 (1996) 6497. 
[18] S. V. Zenkin, Phys. Lett, 395B (1997) 283. 
[19] R. Narayanan, P. Vranas, |hep-lat/ 9702005 . 



6 



[20] P. de Forcrand, M. Garcia Perez, I-O. Stamatescu, |hep-lat/97UT0T2 
[21] P. M. Vranas, |hep-lat/ 9705023 . 



7 



